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OO ' Abstract. We have made use of archival HST BVIJH photometry to constrain the nature of the three discrete 

^ ! sources, Al, A2 and Bl, identified in the double nucleus of NGC 6240. STARBURST99 models have been fitted 

to the observed colours, under the assumption, first, that these sources can be treated as star clusters (i.e. 
single, instantaneous episodes of star formation), and subsequently as star-forming regions (i.e. characterised by 
continuous star formation). For both scenarios, we estimate ages as young as 4 million years, integrated masses 
rangmg between 7 x lO'' Mq (Bl) and 10^ Mq (Al) and a rate of 1 supernova per year, which, together with the 
stellar winds, sustains a galactic wind of 44 Mq yr~^. In the case of continuous star formation, a star- formation 
QQ . rate has been derived for Al as high as 270 Mq yr~^, similar to what is observed for warm Ultraluminous Infrared 

' Galaxies (ULIRGs) with a double nucleus. The Al source is characterised by a mass density of about 1200 Mq 

fT^ , pc~^ which resembles the CO molecular mass density measured in cold ULIRGs and the stellar density determined 

■ in "elliptical core" galaxies. This, together with the recent discovery of a supermassive binary black hole in the 

' double nucleus of NGC 6240, might indicate that the ongoing merger could shape the galaxy into a core elliptical. 

I ' Key words. Galaxies:evolution, Galaxies:individual: interactions, starburst, star clusters 

o ■ 

1. Introduction multiplicity, ULIRGs show circumnuclear bright knots 

L* ^ n , usually identified with concentrations of stars younger 

The importance of starbursts associated with major merg- ^^^^ ^^7. years (e.g., Zepf & Ashman 1993, ScoviUe et 



> 



^ ; ers between galaxies became clear when IRAS revealed 2OOO, Surace et al. 2000, Farrah et al. 2001). At optical 

l2 populations of luminous and ultraluminous infrared galax- , ^i, ttt ron u i • ju 1 

. • /TXT^^ . XTTTT,^ N • , 1 • / p wavclcngths, ULIRGs are characterised by a nuclear spec- 

les (LIRGs and ULIRGs) m the nearby universe (ci. , i,- i, • -ii, c r . n • i j. -^.u j.v, 

■ - - ' ^ ' ^ trum which is either beytert-iike consistent with the pres- 

Neugebauer et al. 1984, Soifer et al. 1984). Interest there- r a/^at r\c<r\ ttatodi-i tjtt • n 

*= ' 1 1 , ence of an AGN or QSO, LINER-hke or Hll-region-like. 

fore remains high in understanding how these starbursts t 4.1, i i j. j. 1 i j / tjtt 

° XTT^^ 1 -11 .In these last two cases, star clusters and/or Mil regions 

are structured. While LIRGs are characterised by an m- ^ u j • ^ r ■ ■ ^- i^j ■^^ j. 

seem to be the dominant source ot ionization (Veilieux et 

frared (IR; 5 — 500^m) luminosity Lir, > 10 L©, ULIRGs , -|qqq\ 

exceed Lir lO^^ L© (Sanders & Mirabel 1996). ULIRGs ^ ' 

contain 10^ - lO^-^ M© of molecular gas (Evans et al. 2001) ^GC 6240 is one of the closest member of the class of 

and are typically the merger of two or more disk galax- ULIRGs, being at a distance of ~ 98 Mpc (for = 75 

ies of nearly equal mass (< 3:1, Sanders 2001). Indeed, s Mpc ). Its far-IR luminosity is (3 - 11) x 10 

NGC 4038/39 (the Antennae), Arp 220 and NGC 7252 (Wright et al. 1984) and it places NGC 6240 at the 

are not only classical examples of LIRGs and ULIRGs, f^int (or "warm") hmit of the ULIRGs class. ^From the 

but also are key benchmarks in the Toomre sequence of P°™t of view of morphology, NGC 6240 is classified as a 

galaxy mergers (Toomre 1977) merger of two massive disk galaxies (Zwicky et al. 1961, 

HST/WFPC2 and NICMOS imaging has revealed that ^osbury & Wall 1979, Tacconi et al. 1999); the colours of 

5% - 20% of the known ULIRGs are multiple mergers, '^^^l^^'^ ^^^^ circumnuclear stellar populations indicate 

probably descending from compact groups of galaxies (Cui '"^^'g^^ possibly been ongoing for the past ^ 

et al. 2001, Bushouse ct al. 2002). Independently of their ^ (Genzel et al. 1998, Tecza et al. 2000). The gravita- 

tional interactions as well as disturbances induced by star 

Send offprint requests to: pasquali@phys.ethz.ch formation are responsible for large-scale dust lanes, loops, 
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shells and tails extending out to ~ 30 kpc (cf. Pasquali 
et al. 2003, hereafter Paper I). The central region of the 
merger is characterised by a double nucleus, with the two 
optically visible nuclei separated by l."5 — 2" on average 
with their apparent separation increasing at shorter wave- 
lengths (Bryant & Scovillc 1999, Scovillc ct al. 2000). The 
H2 and CO emission is found to peak in between the nu- 
clei, possibly in a thick disk structure, and slightly closer 
to the Southern nucleus (van dcr Wcrf et al. 1993, Tacconi 
et al. 1999, Ohyama et al. 2000). This off-set is believed 
to be the result of interactions between the molecular gas 
and an expanding shell centered on the southern nucleus 
(Rieke et al. 1985), or to the gravitational settling of gas 
into a nuclear disk as the interaction evolves (Tacconi et 
al. 1999). 

A remarkable feature of the NGC 6240 nuclei is the 
three discrete sources, Al, A2 and Bl located in the 
southern "A" and the northern "B" nucleus, respectively. 
They were first resolved in HST/FOC B-band images by 
Barbieri et al. (1993) who measured an enhanced blue 
continuum emission corresponding to these sources. In a 
subsequent analysis based on [0II]/[0III]/H/3 line ratios, 
Rafanelli et al. (1997) suggested that Bl is a LINER, and 
A2 a LINER or a HII region (cf. Veillcux et al. 1995). The 
authors also pointed out that the line ratios derived for 
Al are typical of a Seyfert galaxy or a high-excitation HII 
region. A Scyfcrt-like nature would imply the presence of 
an AGN in the A core. These nuclear sources are believed 
to drive a large-scale outflow (Heckman et al. 1987). 

K-band imaging of NGC 6240 reveals the presence of 
red supergiants in the double nucleus (Tccza et al. 2000) 
from which an age between 15 and 25 Myr is inferred for 
the galaxy nuclei. In the radio, 70% of the flux at 20 cm 
is emitted by A and B, of which half comes from two 
unresolved sources and half is diffuse emission (Colbert et 
al. 1994). These compact sources do not align up with Bl 
and A1+A2, but arc shifted westward and less separated, 
probably because of obscuration effects. Their spectrum 
is quite flat {a = 0.6) suggesting that their radio emission 
is powered by the local star burst. 

Beswick et al. (2001) observed neutral hydrogen in 
absorption against the two nuclei of NGC 6240 using 
MERLIN with an angular resolution of 0.2 arcsec. The 
nuclei appear as slightly extended sources at a frequency 
of 5 GHz with a resolution of 0.102x0.055 arcsec and are 
separated by 1.52 arcsec. The radio continuum emission is 
attributed to a combination of AGNs and starburst activ- 
ity. The HI absorption appears to be associated with the 
gas disk that is also seen in CO. Hagiwara, Diamond, & 
Miyoshi (2003) present radio maps of the H2O masers in 
NGC 6240, and show that the masers are concentrated in 
the southern Al nucleus, which evidently contains dense 
gas. 

ISO data have revealed 01 V emission from the nuclei 
which is attributed to an AGN (Genzel et al. 1998, Lutz et 
al. 2003). AGN- like hard X-ray emission is observed below 
10 keV (Iwasawa & Comastri 1998, Vignati et al. 1999). 
High-resolution follow-up observations with Chandra and 



Table 1. The log of NGC 6240 images taken with 
HST/WFPC2 and NICMOS. 



Dataset 


Filter 


Exp. Time in s 


u4ge010.. 


F450W - broad B 


3 X 700.00 


u4ge010.. 


F547M - medium V 


400.00 -1- 800.00 


u4ge010.. 


F814W - I 


3 X 400.00 


n48h09w.. 


FllOW - J 


4 X 39.95 


n48h09... 


F160W - H 


4 X 47.95 


n48h09x.. 


F222M - K short 


4 X 55.94 



the HRC camera in the 0.08 — 10 keV band have been per- 
formed by Lira et al. (2002). In the high resolution X-ray 
map, the emission peaks coincide with Bl and Al. Lira 
et al. used the X-ray fluxes together with optical-to-radio 
luminosities measured for the galaxy's double core and fit- 
ted the resulting spectral energy distribution (SED) with a 
combination of starburst, QSO and blackbody continuum 
spectra. It then turns out that the starburst component 
dominates the observed SED at most wavelengths and the 
dust associated with the starburst is responsible for the 
K-band emission. The AGN accounts mainly for the hard 
X-ray emission and is completely absorbed in the near IR, 
optical and soft X-ray bands. Very recently, Komossa et al. 
(2003) imaged the A and B nuclei with Chandra ACIS-S 
and clearly detected an AGN in each nucleus, so that the 
existence of a binary system of supermassive black holes 
in the core of NGC 6240 has now been firmly established. 

In this paper, we aim to estimate the age, mass and 
reddening of the stellar content of the nuclear sources in 
NGC 6240 and compare these with their properties at ra- 
dio and X-ray wavelengths, in order to: i) discuss the star 
formation mode ongoing in the core of such a gas-rich and 
dynamically extreme system, and ii) possibly shed some 
light on the nature of the end-product of the merging pro- 
cess at work in NGC 6240. 

2. Data processing 

2.1. Observations 

NGC 6240 was observed with HST at several wavelengths. 
Here, we have collected the data acquired with WFPC2 
as part of the GO proposal 6430 (PI van der Marel, cf. 
Gerssen et al. 2001) and NICM0S/NIC2 for the GO pro- 
posals 7219 (PI Scoville) and 7882 (PI van der Werf). The 
image rootnames with their corresponding filters and ex- 
posure times are listed in Table 1. 

We retrieved the WFPC2 and NICMOS images al- 
ready pipeline-processed, i.e. corrected for bias, dark cur- 
rent and flat fleld. We registered the WFPC2 frames 
to the same spatial grid of the images taken with the 
F450W filter, by simply measuring the position of several 
point sources in common and their relative shifts in (X;Y) 
among the available datasets. The (X;Y) shifts were then 
apphed with the IRAF routine IMSHIFT. Once ahgned, 
the images acquired with the same filter were combined 
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Fig. 1. The nuclear sources in NGC 6240, Al, A2 and Bl, 
as they appear in the PC I image (top) and in the NIC2 
H image (bottom). The latter has been scaled to the PC 
angular resolution. 



with the STSDAS task CRREJ to clean them from cos- 
mic rays and have them median-combined to improve their 
S/N ratio. 

The NICMOS pipeline-reduced images were aligned to the 
frames acquired with the FllOW filter with the same tech- 
nique as above. Since they were obtained with Camera 2 
(FOV of 22" X 22" at 0.075"/pix) they essentially over- 
lap with the central area of the PC images containing the 
galactic double nucleus. 

2.2. Multiband photometry 

The high spatial resolution of the Planetary Camera has 
made possible to detect three discrete sources in the nuclei 
of NGC 6240 which were labelled A and B by Barbieri et 
al. (1993). On the basis of their coordinates, our sources 
have counterparts in Barbieri et al.'s FOC images, namely 
Al, A2 and Bl, and in the Chandra images at 0.08 — 10 
keV obtained by Lira et al. (2002). Figure 1 shows the 
F814W and F160W images of these nuclear sources, where 
the NIC2 image in H has been scaled to the resolution of 
the PC I image. 

Their positions were measured with the IRAF routine 
IMEXAMINE, and were used as the input coordinates for 
PHOT in DAOPHOT. Photometry of the nuclear sources 
was performed assuming an aperture radius of 3 pixels 
in the PC images. A median, local background of diffuse 



emission underlying the nuclear sources was computed 
over an annulus of 4-to-6 pixels. 

Aperture fluxes were corrected for the camera's charge 
transfer (in)e0iciency (CTE) following Whitmore & Heyer 
(1998). Aperture corrections to the standard 0".5 aper- 
ture were determined from two bright and point-like iso- 
lated objects in the PC, and subsequently applied to 
the measured fluxes. Only at this stage, fluxes were first 
transformed into the WFPC2 synthetic magnitude sys- 
tem and subsequently into the Johnson magnitude sys- 
tem using the zero points and the colour equations listed 
in Holtzman et al. (1995). 

Because it becomes very difficult to clearly identify 
the Al, A2 and Bl sources in the K band, where emission 
from dust dominates the stellar emission from the nuclei 
(cf. Gerssen et al. 2003), we performed aperture photom- 
etry only on the FllOW and F160W images. We set the 
aperture radius to 2 pixels and measured a median sky 
in an annulus of 5-to-8 pixels. Subsequently, background- 
corrected aperture fluxes were scaled flrst to those corre- 
sponding to the fluxes in a 0".5 radius aperture and then 
to the flux of an infinite radius aperture following the pre- 
scription given in the NICMOS Photometry Cookbook (cf. 
http: /www. stsci.edu /hst /nicmos / performance/ 
photometry / cookbook.html) . 

The BVIJH magnitudes of the nuclear sources and 
their uncertainties are listed in Table 2. 

3. The stellar populations in the nuclear sources 

The Al, A2 and Bl sources have a mean FWHM of 5.2, 
9.6 and 4.2 pixels respectively, a factor of 3 to 6 larger than 
the PC stellar PSF, therefore they are resolved. Prom their 
FWHM we can derive estimates for their tidal radii, which 
vary between 47 and 108 pc. Such values have also been es- 
timated for some of the younger and intermediate-age star 
clusters in the Antennae (Whitmore et al. 1999): knot S 
and cluster #430 have a tidal radius of 450 pc and > 73 
pc, respectively, while the intermediate-age cluster #225 
has a tidal radius of 50 pc. In comparison, all the globular 
clusters in M31 are characterized by Rt < 100 pc (Cohen 
& Freeman 1991, Grillmair et al. 1996), and only two 
Galactic globular clusters have Rt > 200 pc (Djorgovski 
1993). In addition, young and massive star clusters are 
quite often found in the nuclei of late-type spiral and 
barred spiral galax;ies (cf. Sect. 5). On the other hand, 
the nuclear sources in NGC 6240 possibly show spectra 
typical of HII regions (Veilleux et al. 1995, Rafanelli et 
al. 1997) and may well be extended clumps of star forma- 
tion. This option would also be more consistent with the 
findings of an underlying, ~ 1 Gyr-old stellar population 
(Genzel et al. 1998, Tecza et al. 2000), so that the light of 
the newly born OB stars is now dominating the colours of 
the nuclear sources. 

We have fitted the colours observed for the nuclear 
sources [(B-V), (V-I), (V-J) and (V-H) corrected for the 
Galactic E(B-V) = 0.076 mag (Schlegel et al. 1998)] to 
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Table 2. The observed magnitudes of the discrete sources in NGC 6240 nuclei. 



Cluster 


B 


V 


I 


J 


H 


Al 


23.27 ± 0.05 


20.82 ± 0.04 


18.73 ± 0.02 


17.49 ± 0.02 


16.96 ± 0.02 


A2 


21.59 ± 0.03 


19.77 ± 0.04 


18.56 ± 0.05 


17.73 ± 0.03 


17.33 ± 0.03 


Bl 


21.19 ± 0.01 


19.84 ± 0.01 


19.24 ± 0.03 


18.66 ± 0.04 


18.49 ± 0.04 



STARBURST99 (Leitherer et al. 1999) tracks, which were 
computed: 

i) for instantaneous star formation in a star cluster with 
a mass of 10^ Mq and solar metallicity (cf. Sect. 3.1); 

ii) for continuous star formation with a rate (SFR) of 10 
Mq yr~^ and solar metalhcity (cf. Sect. 3.2). 

Indeed, colours depend mainly on reddening and age once 
the metallicity is fixed. In massive clusters stochastic ef- 
fects in the stellar population arc niinmized and so the 
colors become insensitive to the total mass so long as 
the stellar mass function is fixed (Langon & Mouchine 
2000, Cerviho et al. 2002). No metallicity has so far been 
measured for NGC 6240, but McCall et al. (1985) pre- 
dicted a metallicity of about 20% solar. The computed 
STARBURST99 tracks also assume a Salpeter IMF from 
0.1 Mq to 100 Mq, a (Type II) supernova cut-off mass 
of 8 Mq, standard stellar mass loss and theoretical wind 
models. 

The theoretical colours at each time mesh point were 
reddened according to Calzetti's (2001) extinction law, 
and a Xtot value was computed as the sum of the differ- 
ences between the reddened theoretical and the observed 
colours, weighted by the observational uncertainties. For 
each time mesh-point, the reddening was varied from E(B- 
V) = 0.0 to 5.0 mag with a step of 0.01 mag and a set of 
reddening, absolute V magnitude and colours have been 
saved corresponding to the minimum Xtot derived from 
the fit to the observed colours. The lowest value among 
the Xtot minima (Xmin) obtained along the age sequence 
of the evolutionary track has been selected together with 
its relative set of best-fitting reddening, colours and ab- 
solute magnitudes. The latter represent our best-fitting 
model SED to each of the three discrete sources. In addi- 
tion, the reddened theoretical apparent V magnitude was 
scaled to the observed one and the initial track mass or 
SFR was multiplied by this scaling factor in order to es- 
timate either the integrated mass or the total SFR of the 
discrete sources. 

In order to assign a confidence interval to age, mass 
and reddening, we have extracted from the previous fits 
all the (reddening, age, absolute magnitudes, mass) solu- 
tions associated with a Xtot < 1-5 times the Xmin value. 
These are selected to obtain the range of acceptable sub- 
optimal fits to the photometry. We have used these in- 
tervals to construct the probability functions of age, mass 
and reddening for Al, A2 and Bl in the same way as done 
in Paper I for the clusters in the main body and tails of 
NGC 6240. In summary, we combined the ranges spanned 
by the fits to age, reddening and mass into histograms, 
where the value of each bin is the number of acceptable 



solutions. Each bin was then normalized by the total num- 
ber of model solutions for the nuclei of NGC 6240, so that 
a probability function is obtained for the age, mass and 
reddening of the sources. A word of caution applies here, 
whereby a photometry measured through larger apertures 
might give somewhat higher masses and SFRs and might 
negligibly affect ages. 

3.1. The cluster hypothesis 

The best-fitting SEDs obtained under the assumption that 
the nuclear sources in NGC 6240 are star clusters arc 
plotted, as solid lines, in the left-hand panels of Fig. 2. 
Superposed are the sources observed, apparent magni- 
tudes, represented by filled dots. The observational un- 
certainties (cf. Table 2) are smaller than the symbol sizes. 
The x^ obtained for each colour for minimum Xtot is plot- 
ted in the right-hand panels of Fig. 2 with the aim of 
showing which colours agree better with the evolutionary 
tracks and which are discrepant. The lowest Xtot (— 6) 
occurs for the Al source, whose theoretical SED agrees 
with the observations to within the photometric errors 
(except for the J band). The Bl and A2 best-fitting SEDs 
deviate from the observed B band by 0.3 and 0.5 mag 
respectively, and this discrepancy could in part be due 
to contamination by nebular emission (see, e.g., Anders 
& Fritze-von Alvensleben 2003 for possible contaminants 
and their effects on broad-band fluxes as a function of 
age of the dominant stellar population). We have esti- 
mated a line-contribution of ^ 0.1 mag in the F450W 
filter from the integrated central spectrum published by 
Fosbury & Wall (1979) and ~ 0.07 - 0.05 mag in the 
FllOW and F160W filters from the spectra in Simpson 
et al. (1996). Uncertainties in the colour equations used 
to transform WFPC2 magnitudes to the Johnson system 
may also contribute. The best-fitting age, reddening and 
mass of Al, A2 and Bl are listed in Table 3. Adopting 
Calzetti's (2001) extinction law, we estimate for the dou- 
ble nucleus of NGC 6240 an Ay extinction between 4 and 
9 mag, in agreement with the value 3 < Ay < 5 ini- 
tially derived by Fosbury & Wall (1979) and DePoy et 
al. (1986). The estimate given by Lutz et al. (2003) is 
larger. Ay ~ 10 — 15 mag, which has been derived from 
the flux ratio of the Br7 to the [Nell] 12.81/zm lines. The 
authors themselves consider their estimated Ay very un- 
certain because of the difficulty to measure the Br7 line 
of the underlying gas against the stellar component. 

The probability functions (for the sources' age, mass 
and reddening) are plotted in Fig. 3 as dashed histograms, 
together with those derived in Paper I for the clusters 
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Fig. 2. The best-fitting spectral energy distributions (solid lines) obtained for the discrete sources Al, A2 and Bl in 
the nuclei of NGC 6240 with superimposed the observed, apparent magnitudes. The assumption here is that Al, A2 
and Bl are star clusters. In the right panels, the value is plotted for each source and colour for minimum Xtot- 



Table 3. The best-fitting age, reddening and mass for Al, 
A2 and Bl computed under the assumption that these 
sources are stellar clusters. 



Source 


Age (Myr) 


E(B-V) (mag) 


Mass (Mq 


) 


Al 


4.4 ± 0.05 


2.3 ± 0.05 


1.2 (± 0.2) X 


10'' 


A2 


3.9 ± 0.10 


1.7 ± 0.02 


4.2 (± 0.5) X 


10* 


Bl 


4.0 ± 0.05 


1.0 ± 0.02 


2.9 (± 0.2) X 


10^ 



in the main body and tails of NGC 6240. A caveat ap- 
plies here, as discussed extensively in Paper I: the his- 
tograms in Fig. 3 simply give the most probable value(s) 
for the clusters' age, mass and reddening, but do not rep- 
resent the true cluster parameters. Indeed, when the mass 
ranges spanned by the clusters are combined in a proba- 
bility function, their width (i.e. the error bar on the mass 



value most consistent with the data for each cluster) has 
the effect of smearing any intrinsic relationship between 
apparent luminosity and cluster mass, and the probabil- 
ity function may show a peak. In addition to this, evolu- 
tionary effects, such as those associated with a non-coeval 
population of star clusters, also broadens power law Imni- 
nosity functions (cf. Meurer 1995, Fritze-von Alvensleben 
1998, 1999, do Grijs et al. 2001, 2003a). Moreover, lumi- 
nosity selection effects (i.e. limiting magnitude of ~ 25 
mag) prevent the detection of clusters at any age that are 
less massive than 4 x 10^ Mq and allow only ever more 
massive clusters to be revealed at gradually older ages. 
Therefore, the mass distribution of the NGC 6240 clus- 
ters (middle panel in Fig. 3) is heavily compromised by 
the limiting magnitudes for masses lower than 4 x 10^ 
M0, even among the very young star clusters. 
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Because the fits for the nuclear sources are based on 
four independent colours, instead of two as for the other 
NGC 6240 star clusters (cf. Paper I) , the confidence inter- 
vals for Al, A2 and Bl are well determined and narrow: 
they do not overlap, except in age. The discrete sources 
in the nuclei of NGC 6240 certainly stand out because of 
their young ages with respect to the overall cluster popu- 
lation, but, while Bl and A2 fall on the tail of the cluster 
mass and reddening distributions, Al shows extreme val- 
ues for both of these parameters. It is interesting to note 
the increase in E(B-V) from Bl in the northern nucleus to 
Al and A2 in the southern one which is consistent with the 
3D geometry of the double nucleus as derived by Ohyama 
et al. (2003). 

3.2. The hypothesis of continuous star formation 

Only a very small number of compact coeval star-forming 
events with masses like those inferred for the nuclei of 
NGC 6240 in the previous Section have been observed 
(e.g., NGC 7252 cluster W3: Schweizer & Seitzer 1998, 
Maraston et al. 2001; NGC 6745: de Grijs et al. 2003b); 
spectroscopic mass confirmation has only been obtained 
for cluster W3 in NGC 7252 to date. The presence of an 
AGN within both nuclei also points to complex evolution- 
ary histories where we expect circumnuclcar stars to span 
a range in age. We further hypothesize that the circumnu- 
clear SFR is likely to have increased during the ongoing 
merger in NGC 6240. Therefore we have to consider a 
continuous star formation scenario, such as described by 
STARBURST99 models with a SFR of 10 Mq yr"! and 
solar metallicity. The resulting best-fitting model SED is 
shown in Fig. 4 for each source, where we also indicate 
the associated with the observed colours for minimum 
Xtot> ^^'i the best-fitting properties of the three nuclear 
regions are listed in Table 4. Similarly as for the cluster 
scenario, we determine for the galaxy double core an Ay 
extinction ranging from 3 to 9 mag. 

The minimum Xiot particularly large for Bl and A2 
as it is dominated by the large discrepancy in the (B-V) 
colour (0.7 and 0.5 mag, respectively). As already sug- 
gested, this deviation could in part be produced by con- 
tamination from nebular emission in the B band; uncer- 
tainties in the colour equations used to transform WFPC2 
magnitudes to the Johnson system may also contribute. In 
the case of Bl and A2, the minimum Xtot is a factor of 
2 worse than obtained under the assumption that these 
sources are star clusters. The SEDs derived for Al for the 
two scenarios are similar (xfot of 5.91 for a cluster SED 
against 6.56 for a continuous star-forming region). The 
difference in (B-V) and (V-I) between the evolutionary 
tracks employed so far (with a fixed cluster mass of 10^ 
Mq and for a SFR of 10 Mq yr~^) are within the observa- 
tional errors for ages of few million years. Therefore, given 
the observed colours of Al, A2 and Bl, we are not able to 
discriminate between a cluster and a sharp SBST model. 
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Fig. 3. The probability distributions of age and mass in 
logarithmic units, and reddening on a linear scale. The 
solid line represents the probability distribution for the 
population of star clusters detected in the NGC 6240 main 
body, the dashed line traces the clusters identified in the 
galaxy's tails (cf. Paper I). The dashed bins show the most 
probable values from fits to the Al, A2 and Bl sources. 

Probability functions have also also derived and are 
plotted as dashed histograms in Fig. 5, together with the 
probability functions computed for the NGC 6240 clusters 
(from Paper I). 

Once again, the discrete sources in the nuclei of 
NGC 6240 are among the youngest stellar systems in the 
galaxy, with ages between 3.5 to 5.4 x 10^ years, thus 
implying that young stars dominate the circumnuclear re- 
gions in optical. The integrated mass of Al, A2 and Bl 
overlaps with the high-mass tail of the overall NGC 6240 
cluster distribution, with Al showing an extreme value 
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Table 4. The best-fitting age, reddening, SFR and integrated mass for Al, A2 and Bl computed under the assumption 
that these sources are star-forming regions and governed by continuous star formation. 



Source 


Age (Myr) 


E(B-V) (mag) 


SFR (Mq yr-^) 


Integrated Mass (M©) 


Al 


4.4 ± 0.9 


2.4 ± 0.07 


266 ± 109 


1.3 (± 1.0) X 10^ 


A2 


4.3 ± 0.4 


1.5 ± 0.05 


21.5 ± 5.7 


9.6 (± 3.4) X lO'' 


Bl 


4.3 ± 0.7 


0.8 ± 0.08 


1.5 ± 0.6 


7.2 (± 3.7) X 10^ 




Fig. 4. As Fig. 2, but now assuming a continuous star formation scenario. 



once again. This same trend can be recognized in the red- 
dening histogram. 

3.3. Dependences on metallicity 

To test the robustness of our results against metallicity, we 
have repeated the above analysis using STARBURST99 
tracks computed under the same assumptions, but with 
the metalhcity set to Z = 0.008 (LMC-hke). Tables 5 and 
6 summarize the best-fitting parameters obtained for the 



cluster and the continuous star formation scenarios, re- 
spectively. 

We conclude that the properties derived for a solar- and a 
LMC-like metallicity are well within the uncertainties of 
the best fits for each of the scenarios modelled. A bet- 
ter overall agreement between the solar and LMC-like 
metallicity cases is achieved for the cluster hypothesis. In 
the case of continuous star formation, the sources' ages 
and star-formation rates arc systematically larger for Z = 
0.008 than for Z = Zq. Although the Xtot values of the 
best fits are comparable, an LMC-like metallicity intro- 
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Table 5. The best-fitting age, reddening and mass for Al, A2 and Bl computed under the assumption that these 
sources are stellar clusters with an LMC-like metallicity, Z = 0.008. 



Source 


Age (Myr) 


E(B-V) (mag) 


Mass (Mq) 


Xtot 


Al 


4.7 ± 2.0 


2.4 ± 0.2 


2.6 (± 1.4) X 10^ 


6.39 


A2 


4.3 ± 1.5 


1.6 ± 0.2 


4.2 (± 2.8) X 10® 


18.87 


Bl 


4.6 ± 1.6 


0.9 ± 0.2 


2.0 (± 1.2) X 10^ 


57.05 



Table 6. The best-fitting age, reddening, SFR and integrated mass for Al, A2 and Bl computed under the assumption 
that these sources are star-forming regions and governed by continuous star formation. Metallicity is now set to Z = 
0.008. 



Source 


Age (Myr) 


E(B-V) (mag) 


SFR (Mo yr-^) 


Integrated Mass (Mq) 


Xtot 


Al 


6.1 ± 2.5 


2.4 ± 0.08 


412 ± 219 


2.5 (± 1.8) X 10^ 


5.86 


A2 


6.2 ± 2.7 


1.5 ± 0.13 


33.4 ± 21.6 


2.1 (± 1.7) X 10** 


.33.64 


Bi 


(>.:> ± 


0.8 ± 0.21 


2.8 ± 2.:i 


i.8 (± 1.6) X iO" 


8.'). 73 



duces larger uncertainties on the parameters than a solar 
chemical composition. Hereafter, we will adopt the results 
of the evolutionary tracks with Z = Zq. 

4. Discussion 

4.1. Nuclear stellar populations 

A starbust can be modeled either as a single and instan- 
taneous event of star formation, or as an act of continuous 
star formation extended over a short period of time. We 
have computed a STARBURST99 model for both scenar- 
ios, and assuming a total cluster mass of 10^ Mq in the 
case of a single stellar burst and a SFR of 10 Mq yr~^ 
for the continuous star formation hypothesis. We have fit- 
ted the photometry measured for the multiple nuclei of 
NGC 6240 Al, A2 and Bl to the synthetizcd populations, 
hoping that the observed colours could discriminate be- 
tween the two scenarios. Within the observational uncer- 
tainties, the results almost overlap and do not support 
either of the two models more strongly than the other. 

4.1.1. Ages 

Both models indicate that the optical SEDs of Al, A2 and 
Bl are dominated by light from stellar populations with 
ages of ~ 4 Myr and moderate-to-heavy obscuration. As 
already discussed by Lira et al. (2002), these optical SEDs 
arc not uniquely tied to the presence of the AGN in either 
Bl or Al. 

We notice that all three nuclear condensations in 
NGG 6240 resemble other Star burst/ AGN systems, where 
large contributions to the light come from young stel- 
lar populations. Indeed, Gonzalez Delgado et al. (2001) 
found that young and intermediate-age stars are present 
in the nuclear and circumnuclear regions of 45% of Seyfert 
2 galaxies. The contribution of young and intermediate- 
age stars to the nuclear stellar population in Seyfert 2 
galaxies can be as high as 50% and 40% respectively (e.g. 
NGC 5135 and NGC 7130; Gonzalez Delgado et al. 2001). 



4.1.2. Stellar Mass Estimates and Star Formation 
Rates 

Under the assumption of a single and instantaneous burst 
of star formation, we have derived stellar masses of about 
1 X 10^, 4 X 10^ M0 and 3 x 10^ M© for Al, A2 and 
Bl, respectively. If the cluster hypothesis were correct, 
these nuclei would be the most massive known, exceed- 
ing by a factor 5 — 10 the most massive globular clus- 
ters. In particular, Bl would be the counterpart of the 
W3 cluster in NGC 7252 (at a distance of ~ 15" from the 
galactic centre; Whitmore et al. 1993) for which Maraston 
et al. (2001) estimate a total mass of 3.7 x lO"^ Mq 
and an age of about 300 Myr, using high-resolution spec- 
troscopy. Another massive cluster has been detected in 
NGC 1316 with a mass of about 1.4 x lO"^ M© and an age 
of 3 Gyr (Goudfrooij et al. 2001). In addition, dc Grijs 
et al. (2003b) have extended the cluster mass function in 
the interacting galaxy NGC 6745 up to masses of several 
xlO* Mq, although their mass estimates are thus far only 
based on multi-passband HST photometry. The masses of 
Al, A2 and Bl, combined with the clusters mass distri- 
bution in the outer (tails) and inner regions of the galaxy 
(cf. Paper I), would point to a radial gradient in the clus- 
ter mass, whereby less massive clusters would preferen- 
tially form far from the galactic centre and more mas- 
sive clusters closer. In the assumption that the dynamical 
timescale of the merger is larger than the age of the nuclear 
sources, we would expect to find more massive clusters in 
the galactic nuclei not as the result of mass segregation, 
but as a consequence of larger amounts of gas funneled to- 
wards the galactic centre by the merger dynamics. Indeed, 
the central location of such massive clusters as Al, A2 
and Bl is particularly interesting. So far, there have been 
a number of detections of a single massive cluster in the 
centre of a galaxy. For example, Carollo et al. (1997) re- 
solved compact sources in the nuclei of 18 spiral galaxies 
(with an exponential bulge) in a total sample of 35 galax- 
ies, for which a mass of 10^ — 10^ Mq and an age > 1 Gyr 
are inferred (Carollo et al. 2001). More recently, Boker et 
al. (2002) identified a compact cluster at the centre of 59 
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Fig. 5. As Fig. 3, but under the assumption that Al, A2 
and Bl are star-forming regions. 

late-type spiral galaxies in a sample of 77 galaxies. For 
example, the cluster at the centre of NGC 4449 was dated 
between 6 and 10 Myr old and as massive as 4 x 10^ 
(Boker et al. 2001). Colina et al. (2002) identified a cluster 
of 4 Myr old and as massive as 10^ Mq in the nucleus of 
NGC 4303, which is a barred spiral and thought to be a 
low-luminosity AGN. By analogy to what is observed for 
NGC 6240 (cf. Lira et al. 2002), this nuclear cluster dom- 
inates the optical emission of the galaxy core. Additional 
clusters, although slightly yoimger (3 — 3.5 Myr) and less 
massive (7 — 8 x 10^ M©), have been detected in the cir- 
cumnuclear regions of NGC 4303 (Colina et al. 2002). 

Assuming continuous star formation, we obtained for 
Al, A2 and Bl a SFR of 266, 21.5 and 1.5 Mq yr"!, re- 
spectively, which gives an integrated mass of about 1 x 
10^, 1 X 10* and 7 x 10^ M© respectively. Star formation 
appears to be by far more active in the southern nucleus, 
A, which is consistent with the southern nucleus being 
brighter at radio and X-ray wavelengths, in H2, CO and 
H2O maser (Hagiwara et al. 2003) emission. This might 
point to a larger reservoir of gas in the southern nucleus 



compared to the northern on (see Sect. 4.2 for a discussion 
on the relation between the star formation history and gas 
supply in NGC 6240). In particular, the SFR estimated for 
Al in the optical is of the same order as the SFR derived 
from the integrated IR luminosity of NGC 6240 (~ 140 
Mq yr~^, Hcckman ct al. 1990), thus indicating that the 
nuclear sources dominate the galaxy's IR emission and 
star formation activity. For these integrated masses and 
for the cluster masses derived above, we estimate a total 
bolometric luminosity of ~ 10^^ L©. Lutz et al. (2003) 
have computed Lboi — 5 x 10^^ L© from comparison with 
starburst template spectra in the [Nelll] 15.55/im / [Nell] 
12.82/Km, with the caveat that this may be an underesti- 
mate since NGC 6240 is more obscured than the template 
star bursts. 

The SFR in the double nucleus of NGC 6240 has also 
been measured at radio wavelengths by Beswick et al. 
(2001). They have derived, from the 1.4 GHz luminosity, 
a SFR of about 83 M© yr'^ adopting a Salpeter IMF 
and assuming that stars initially more massive than 8 
Mq become supernovae (SNe) with significant radio emis- 
sion. This value is a factor of 2 to 5 lower than what 
we have estimated for Al, using BVIJ photometry and 
STARBURST99 models based on a Salpeter IMF with 
only stars more massive than 8 Mq exploding as type 
II SNe. This discrepancy could have a variety of possi- 
ble sources, including the calibration of the techniques, 
possible impact of rapid evolution in NGC 6240, and the 
offset between regions that are bright in the optical and 
the radio, as evidenced by the radio centres being shifted 
westward and having a smaller projected separation than 
the optical. This smaller distance optical light and radio 
wavelengths, with the radio centres being shifted west- 
ward and less separated. This smaller distance is probably 
an effect of extinction increasing in between Al and Bl 
(cf. the 3D structure of the double nucleus as modeled by 
Ohyama et al. 2003). Therefore, eventual star-forming re- 
gions between Al and Bl would be hidden from detection 
at optical wavelengths and become visible only at radio 
wavelengths. As a caveat, the radio SFR could in princi- 
ple be associated with star-forming regions different from 
Al, A2 and Bl, which are not detected at optical wave- 
lengths because of their high extinction. However, despite 
these uncertainties, the point remains that NGC 6240 has 
a prodigious SFR in its double nuclei. 

In terms of SFRs and ages, the nuclear sources in 
NGC 6240 closely resemble those detected in Arp 299. 
The latter appear to be 4 - 7 Myr old, with SFRs be- 
tween 40 and 140 Mq yr~^ and integrated masses in the 
range (0.3 - 7) x 10* Mq (i.e. 1 to 2 orders of magni- 
tude more massive than star clusters and HIT regions in 
Arp 299, Alonso-Herrero et al. 2000). According to Ballo 
et al. (2003), Arp 299 harbours an active AGN in both 
nuclei and would be, with NGC 6240, the only merging 
system with two AGNs. Although the statistics are very 
poor, it would seem that the more powerful AGNs may 
be coupled to larger nuclear SFRs and larger total nu- 



10 



Anna Pasquali, John S. Gallagher and Richard de Grijs: Nuclear star formation in NGC 6240 



clear masses in newly born stars (~ 4 Myr old), as already 
pointed out by Gonzalez Delgado et al. (1998). 



4.1.3. The High Energy Budget 

Irrespective of the scenario adopted (star clusters or con- 
tinuous star formation), the STARBURST99 code com- 
putes, for the inferred ages and total masses of Al, A2 
and Bl a SN rate of about 1 SN yr"-'^, in agreement with 
what Colbert et al. (1994) and Beswick et al. (2001) es- 
timated from the nuclear radio emission at 20 cm and 
1.4 GHz, respectively. The largest contribution to this SN 
rate comes from Al. These SN explosions would power a 
galactic superwind which manifests itself as a gas outflow 
in the Ha and Br7 images (cf. Heckman et al. 1987, Paper 
I)- 

The predicted mechanical energy flux released by the 
stellar winds and SN explosions in Al, A2 and Bl is ~ 6 
X 10^^ ergs s^^ from both the continuous and instanta- 
neous star-formation scenarios. This value is in agreement 
with the predictions of Heckman et al. (1990). Indeed, if 
a total IR luminosity of (3 — 11) x 10^^ L© is assumed 
for NGC 6240, equation (11c) in Heckman et al. gives an 
energy flux from SN explosions and stellar winds in the 
range (2 — 8) x 10^'' ergs s^^. According to the superbub- 
ble model by Mac Low & McCray (1988), a conversion 
factor of total wind/SN energy flux to kinetic energy of 
the gas of ^ 30% implies a kinetic energy of < 2 x 10'*'' 
ergs s~^. Adopting a velocity spread of the Ha filaments 
of 1000 km s~^ (Heckman et al. 1990), we have computed 
a galactic wind, mass-loss rate of about 44 Mq yr~^. 

How do the STARBURST99 models compare with the 
observed, soft X-ray emission in NGC 6240? Schulz et al. 
(1998) determined, from ROSAT observations, a soft X- 
ray luminosity of about 10*^ ergs s^^. Based on the super- 
wind model by Mac Low & McCray (1988) and Heckman 
et al. (1996), Schulz et al. suggested that this luminosity 
is due to an input energy flux of 10*^ ergs s~^, equivalent 
to 3 SNe per year. Such a flux can expand a superbubble 
out to a radius of 10 kpc in 3 x 10*" yr. The calculations 
assumed an ISM density of 0.1 cm^"^ and a filling factor 
of the X-ray emitting gas of about 0.2. We have repeated 
the exercise of Schulz et al. by assuming an input energy 
of 6 X 10^3 ergs s"! (from the STARBURST99 fits to the 
observed photometry of Al, A2 and Bl with Z = Zq) and 
a kinetic energy of 2 x 10^^ ergs s~^. From Schulz et al.'s 
equation (3), we estimate a soft X-ray luminosity of ~ 6 
X lO'^^ ergs s~^, which agrees with the observed X-ray lu- 
minosity within a factor of < 2. The agreement strongly 
depends on the uncertainties in the observations and the 
model assumptions. However, it could be improved if, for 
example, the conversion factor of input energy flux to ki- 
netic energy were higher than 27% (up to perhaps 50%) 
and/or the filling factor of the superbubble were smaller 
than 0.2, by perhaps a factor of 2. 



4.2. Star Formation Lifetimes and Histories 

The nuclear condensations Al, A2 and Bl appear to be 
"embedded" in a circumnuclear population of ~1 Gyr old 
(Genzel et al. 1998, Tccza et al. 2000), possibly implying 
that the NGC 6240 nuclei have been undergoing a large 
spread in the star formation activity. This co-existence of 
old and young stellar populations has also been seen in the 
Antennae nuclei, where Berta et al. (2003) have been able 
to date two star-formation episodes; one occurred about 1 
Gyr ago and the other is still ongoing. The authors suggest 
that the bimodal star-formation history in the Antennae 
nuclei is the result of a double encounter in the dynami- 
cal evolution of the galaxies' merger. Similar star forma- 
tion patterns are seen in other classes of interacting galax- 
ies; e.g. the Seyfert 1 NGC 3227 (Schinnerer et al. 2001). 
However, it is also possible that the pre-merger systems 
could provide a substantial fraction of the older stars seen 
in NGC 6240. The correlation between nuclear black hole 
mass and the stellar mass of the bulges of non-interacting 
galaxies implies that the pre-merger galaxies in NGC 6240 
both had significant bulges (Kormendy & Richstone 1995, 
Marconi & Hunt 2003). 

The data also indicate that substantial variations in 
SFRs occur over very short time spans, < 10 Myr, in 
mergers. For example, in Arp 299 Alonso-Herrero et al. 
(2000) detected an extended period of star formation that 
began > 7 Myr before the most recent episodes with ages 
of ^ 4 Myr. The extremely young stellar population ages 
of ~ 4 Myr that we find for the three NGC 6240 nuclei is 
another indication of this effect. Rapid variations in SFRs 
might be expected to be associated with the compact, high 
density central gas concentrations that are characteristic 
of ULIRGs and related classes of starburst galaxies (e.g. 
Solomon et al. 1997). These circumnuclear regions have 
sizes of 0.1 — 0.5 kpc and internal velocities of > 100 km 
s~^ and associated dynamical scales of only ~ lO'' yr. It 
should be no surprise that star formation near the cen- 
tres of merging galaxies is dynamic, but given the huge 
gas concentrations in these regions, typically ~ 10^ M©, 
the effects of SFR fluctuations are spectacular. Given that 
our mass estimates for Al, A2 and Bl are comparable 
to within a factor of < 2 between the instantaneous and 
continuous star formation scenarios, our results are fairly 
robust with respect to the star formation history adopted. 

The time to produce the stellar mass in a system if the 
current SFR remains constant is an useful indicator of evo- 
lutionary rates. From Table 4 we see that this scale is ^ 
10^ yr for all the three nuclear condensations. Another es- 
timate of the expected time scales for the evolution of the 
NGC 6240 nuclei can, in principle, be obtained by exam- 
ining the mass balance of the nuclear interstellar medium, 

MiSM = Ma - SFR - Mwind, 

where Ma is the gas accretion rate into each nucleus and 
Mwind is the mass loss rate out of each nuclear region 
fueling the galactic wind. Using nucleus Al as an example, 
we have M^md < 40 Mq yr^^ (cf. Sect. 4.1.3 and Heckman 
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et al. 1990) and SFR ~ 270 Mq yv'^. The major unknown 
is the accretion rate into the nucleus. We can estimate this 

rate by A/g^ ~ -A^g,nuclcar/^cross, nuclear whcrC ^cross, nuclear 

R/V ~ 10'' yr. Adopting a nuclear gas mass of 10^ Mq 
from Tacconi et al. (1999), we then estimate Ma ~ 100 
Mq yr-1. 

The characteristic time scale for the starburst in nu- 
cleus Al is then Tburst ~ —^g I < 10^ yr. The combi- 
nation of high SFRs, large flow velocities, and small sizes 
of the relevant regions demand rapid variations in SFRs. 
Even if one invokes a maximal estimate of the total mass of 
interstellar gas in NGC 6240 of 3 x 10^° M© (Gcorgakakis 
et al. 2000, allowing for the mass of **He) and a short dy- 
namical time scale of 100 Myr for the entire system, we 
see that the current starburst can be sustained for < 100 
Myr. The actual lifetime of a burst with the current SFR 
is probably considerably shorter under more realistic as- 
sumptions as to how much gas actually can fuel the nuclear 
starburst region. 

4.3. Evolution of the Nuclei 

4.3.1. Compact Disks and the AGN 

While the two AGN in NGC 6240 and their associated star 
clusters appear to be simple objects at our angular resolu- 
tion, kinematic studies supply a glimpse of their true com- 
plexity. Observations of the central stellar velocity field 
by Tecza et al. (2000) from the ground in the K-band 
show pronounced and organized radial velocity gradients 
around nuclei Al and Bl. They interpret these velocity 
patterns as resulting from the rotation of inclined stellar 
objects, such as disks or flattened bulges, and by apply- 
ing this model derive dynamical masses of >10^ Mq. High 
angular resolution observations of emission line velocities 
and widths by Gerssen et al. (2003) show more complex- 
ity. Their results could arise from some combination of 
nuclear disks and large-scale gas flows associated with the 
starbursts and ongoing interaction. Disturbed large-scale 
flows in the gas that are largely decoupled from the nu- 
clei are also consistent with what is seen in CO emission 
(Tacconi et al. 1999), H2 IR emission lines (Ohyama ct al. 
2000, Ohyama et al. 2003) and in HI absorption against 
the AGNs (Beswick et al. 2001). 

Despite the complex gas kinematics seen towards Al 
and Bl, it is likely that both AGNs are surrounded by 
compact gas disks that have substantial rotational sup- 
port. This model is consistent with the K-band stellar 
kinematics, and also with the presence of OH maser emis- 
sion from Bl (Hagiwara, Diamond, & Miyoshi 2003). 
Nuclear maser emission in ULIRGs is thought to occur 
in rotating disks of dense gas with sizes of a few hun- 
dred pc, which provide the large column densities of gas 
at near constant velocity that are required for masers (cf. 
Philstrom et al. 2001). When other nuclear starbursts with 
AGNs are studied with suflncient precision, molecular gas 
disks with radii of ~100 pc and masses >10^ Mq are found 
(e.g., Mrk 231; Bryant & Scoville 1996, Arp 220; Scoville 



et al. 1997, NGC 4303; Cohna & Arribas 1999). While 
such structures seem to be generic features of strong in- 
teractions involving gas-rich galaxies, it is not clear which 
physical properties determine their size and mass scales 
(e.g., Downes & Solomon 1998, Jogee et al. 2002). 

4.3.2. AGN and Pre-Merger Bulge Masses 

The combined total luminosity of the two AGNs in 
NGC 6240 is estimated by Lira et al. (2002) to be ~5 x 
10^^ erg s"-'^, comparable to what is produced by the star- 
burst (see also Lutz et al. 2003). We therefore assume that 
each nucleus produces ^10^^ L©. A lower boimd estimate 
for the black hole masses of A/bh ^10^"'' Mq then fol- 
lows from the Eddington limit argument. Assuming that 
most of this mass was deposited in the black holes before 
the present merger, we can estimate properties of the pre- 
merger system bulges on the basis of empirical correlations 
between Mbh and host galaxy bulge properties. Thus we 
find that the pre-merger bulge velocity dispersion would 
have been «150 km s~^ and the bulge B-band luminosi- 
ties about 10^'^ Lq corresponding to a stellar mass near 
10^" Mq. The presence of dual AGN in NGC 6240 is con- 
sistent with it being a major merger involving two giant 
spiral galaxies. 

The similarity between conditions in Arp 299 and those 
in NGC 6240 suggests another possible evolutionary trend. 
NGC 6240 and Arp 299 are the only two merging systems 
known to contain dual AGNs. In both cases, it appears 
that the amount of mass going into the circumnuclear 
starbursts considerably exceeds that being accreted by the 
nuclear black holes. Since energy production from accre- 
tion onto a black hole is more than order of magnitude 
more efficient than what is provided by nuclear energy in 
stars, and the AGN and starburst luminosities are com- 
parable in NGC 6240, we predict a black hole accretion 
rate of ^ 3 — 30 M© yr~^ in nucleus Al. From our earlier 
discussion of the lifetime of the starburst, we take 30 Myr 
as a generous estimate of the main starburst lifetime. We 
then see that if the system continues in its present state, 
the black holes could grow in mass by >10^ Mq prior to 
the eventual merger of the binary nuclei into a single ob- 
ject. Our analysis thus suggests the nucleus of the final 
NGC 6240 merger will have Mbh >2x10^ Mq. 

4.3.3. The A2 Cluster 

Our photometric study demonstrates that nucleus A2 is 
also a huge concentration of young stars. However, it 
shows no evidence for being an AGN, and therefore could 
be an example of an extraordinarily luminous star clus- 
ter. For example, young super star clusters with similar 
masses are found in the centre of the Arp 220 merger, 
where the extraordinary N2 and N3 clusters described by 
Shioya, Taniguchi, & Trentham (2001) are not necessarily 
associated with nuclei. The fate of such clusters is some- 
what uncertain. However, due to their high masses, these 
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clusters likely will have short lifetimes driven by the dy- 
namical friction time scale for them to spiral into the cen- 
tres of their host galaxies, where they will merge with the 
AGNs/nuclear star clusters. This type of event can help 
to populate the bulge of the remnant galaxy (cf. Noguchi 
1999, Kim, Morris, & Lee 1999). 

4.3.4. The Fate of NGC 6240 

The Bl, Al and A2 integrated masses give a mass den- 
sity of about 16.0, 1.2 x 10^ and 19.0 Mq pc'^, respec- 
tively, as measured in a circular aperture of radius 0".l, 
and assuming continuous star formation. For a single and 
instantaneous burst of star formation, the cluster mass 
density increases to 67.0, 1.4 x 10^ and 80.0 M© pc~^. In 
either scenario, these values are similar to the mass den- 
sity of the molecular gas (10^ — 10^ M© pc~^) integrated 
over the central R < 0.5 — 1 kpc of most nearby ULIRGs 
(Sanders & Mirabel 1996), where 0.5 kpc translates into 
1" at the distance of NGC 6240, i.e. the entire double core 
of NGC 6240. Evans et al. (2002) estimated the core mass 
density in molecular gas for a sample of ULIRGS with 
double nuclei and obtained values of 970 to 8.6 x 10^ Mq 
pc~^ per nucleus. Compared with this sample, Al is simi- 
lar to IRAS 12112+0305SWand IRAS 14348+1447SW. In 
addition, the mass density of Al falls in the low-mass den- 
sity tail of the distribution of early-type galaxies, whose 
masses were measured in an aperture with a radius of 
0".l by Faber et al. (1997). Specifically, Al appears to fit 
in with the "core" galaxies of Faber et al.'s sample. Faber 
et al. (1997) separate early- type galaxies into two cate- 
gories: the power-law galaxies with a fairly steep surface 
profile in their surface brightness, extending inward to the 
smallest resolvable radius, and the core galaxies showing a 
broken power-law profile which changes to a significantly 
shallower slope at a specific break radius. Core galaxies 
are bright objects with My < -20.5 (H,, = 80 km s^^ 
Mpc^^) and are found in clusters as well as in the field. 

The mass densities derived for Bl, Al and A2 sug- 
gest that, whatever the nature of the burst is, if their fate 
is to merge together, then Bl, Al and A2 would likely 
form the core of an elliptical or at least a large bulge of 
an early-type spiral galaxy, which NGC 6240 could evolve 
into, once the merger of the parent galaxies is completed. 
The findings of Komossa et al. (2003) support this sce- 
nario. The presence of a heavily obscured AGN in the A 
nucleus has long been known, but the definitive detection 
of an AGN in the B nucleus, hence of supermassive binary 
black holes, comes from the very recent Chandra ACIS- 
S data. According to Milosavljevic & Merritt (2001) and 
Milosavljevic et al. (2002), during a galactic merger the 
two black holes fall to the centre of the merger and form 
a bound pair. Their spiraling towards each other releases 
their binding energy to the surrounding stars, modifying 
the central stellar density of the merger. Consequently, 
the initial central cusp in the stellar density evolves into a 
shallower distribution, p ~ r~^, similar to the luminosity 



profile observed at the centre of elliptical core galaxies. 
This evolution occurs within 10^ — 10^ yr after the binary 
black hole has formed; the cusp continues to flatten there- 
after as the binary ejects stars via gravitational kicks away 
from the nucleus of the merger. The process should ter- 
minate when the binary reaches a separation of between 
0.01 and 1 pc. 

The black holes detected in the double nucleus of NGC 
6240 are about 1.5 arcsec apart (Komossa et al. 2003), 
which translates into 700 pc at a distance of 98 Mpc, there- 
fore NGC 6240 may be at the beginning of the "cusp co- 
alescence" defined by Milosavljevic & Merritt (2001) and 
may transform into an elliptical core galaxy. This agrees 
with the conclusion by Genzel et al. (2001) who mea- 
sured a rotational velocity and a velocity dispersion in 
NGC 6240 consistent with the properties of elliptical and 
lenticular systems. 



5. Summary 

We have analysed HST data available for three discrete 
sources identified in the double nucleus of NGC 6240. We 

have fitted the observed colours to STARBURST99 mod- 
els to estimate their ages, masses and intrinsic redden- 
ings. Our fitting technique (cf. Paper I) relies on observed 
colours corrected for the Galactic extinction in the direc- 
tion of NGC 6240 and develops through five steps: i) the 
synthetic colours at each time mesh-point of the model 
are reddened by E(B-V) increasing from 0.0 to 5.0 mag 
in steps of 0.01 mag; ii) at each step in E(B-V), a Xtot is 
computed as the sum of the colour differences between ob- 
servations and evolutionary tracks weighted by the obser- 
vational errors and squared; Hi) for each time mesh-point 
of the model, a minimum Xtot extracted from all the Xtot 
values obtained by varying the reddening and the model 
age, magnitudes and colours; iv) once all the time mesh- 
points of the model have been fitted, the smallest Xtot 
(Xmin) identified among all Xtot minima and its corre- 
sponding E(B-V); model age, magnitudes and colours are 
assigned to the cluster concerned; v) the selected model 
magnitudes are reddened by the final selected E(B-V), cor- 
rected for the distance modulus to NGC 6240 and scaled 
to the observed apparent magnitudes. The assumed SFR 
is multiplied by this magnitude scaling factor to derive the 
actual SFR of the sources. We have made use of two sets 
of STARBURST99 models, computed for two different 
metallicities, solar- and LMC-like. The results obtained 
from these two sets of evolutionary tracks agree reason- 
ably well and are within the uncertainties of our fitting 
procedure. Therefore, the subsequent analysis and discus- 
sion of the properties of the nuclear sources have been 
based on the STARBURST99 models with solar metal- 
licity. The nuclear condensations in NGC 6240 have four 
independent colours available, (B-V), (V-I), (V-J) and (V- 
K), which allow us to determine a unique combination of 
age, reddening and SFR with relatively small uncertain- 
ties. 
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The probability distributions for the nuclear sources 
in NGC 6240 have been derived from the fits with Xtot ^ 
1.5 times the x^in vahic of the fitting procedure. Under 
the assumption that these sources are clusters, they are 
dated as young as 4 Myr and are heavily obscured, with 
E(B-V) = 1 — 2.3 or Ay = 4 — 9 in agreement with spectro- 
scopic measurements. Their mass estimates, however, are 
extremely high, from 3 x lO'^ to 10^ M0. A mass of lO'^ 
Mq, as foimd for Bl, has been so far determined and con- 
firmed spectroscopically only for cluster W3 in NGC 7252, 
while a lower mass of a few 10^ Mq has been estimated 
for the nuclear clusters in NGC 4303 and NGC 4449. On 
the other hand, the tidal radius of the nuclear sources in 
NGC 6240 is larger than measured for bona-fide compact 
clusters and thus this favours the earlier suggestion that 
these sources are giant Hil regions instead. In this case, a 
STARBURST99 model with continuous star formation is 
probably more appropriate, which results in: 

1) young ages of about 4.3 x 10^ yr and reddenings in the 
range 0.8 < E(B-V) < 2.4 mag; 

2) SFRs between 1.5 (Bl) and 266 (Al) M© yr-\ which 
define integrated masses from 7 x 10^ to 10^ M©; 

3) mass densities from 16 (Bl) to 1200 (Al) Mg pc'^. 

4) a SN rate of 1 SN yr'^ and a galactic wind of 44 M© 
yr~^. Comparison with previously determined SN rates 
and galactic wind estimates based on NGC 6240's Lir, 
indicates that the nuclear sources dominate the galaxy's 
total IR luminosity and star formation activity. 

Irrespective of the scenario assumed, the properties of 
the nuclear sources are such to likely sustain the observed 
soft X-ray emission. They can also be used to estimate a 
mass of ~ 10^° M© for the bulge of the parent galaxies and 
a black- hole mass of > 2 x 10*^ M© for the end-product 
of the coalescence of the binary black hole. 

Star formation rates such as those found for Bl and 
A2 are commonly observed in nearby blue compact galax- 
ies such as NGC 7673 and Lyman break galaxies, while 
the SFR value derived for A2 is typical of warm ULIRG 
galaxies with double nuclei. Similarly, the mass density 
in Al is consistent with the mass densities of molecular 
gas in ULIRGs and the stellar mass densities measured 
in elliptical core galaxies. The supermassive binary black 
holes recently discovered in NGC 6240 (Komossa et al. 
2002) could in time transfer energy to the stars in Al, 
A2 and Bl, thus lowering the central stellar density in 
NGC 6240 nuclei to a flat p ~ profile, typical of ellip- 
tical core galaxies (cf. Milosavljevic & Merritt, 2001). The 
balance between the relative efficiencies of star formation, 
dynamical friction, and the slingshot effect of the coalesc- 
ing binary black hole will ultimately determine whether 
the merger in NGC 6240 will produce a cusp or a core 
elliptical. 
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